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We present an analytical calculation of the velocity of a single 180° domain wall in a magnetic 
structure with reduced thickness and/or lateral dimension under the combined action of an exter- 
nal applied magnetic field and an electrical current. As for the case of field-induced domain wall 
propagation in thick films, two motion regimes with difi'erent mobilities are obtained, below and far 
above the so-called Walker field. Additionally, for the case of current induced motion, a Walker-like 
current density threshold can be defined. When the dimensions of the system become comparable 
to the domain wall width, the threshold field and current density, stating the wall's internal struc- 
ture stability, are reduced by the same geometrical demagnetising factor which accounts for the 
confinement. This points out the fact that the velocity dependence over an extended field/current 
range and the knowledge of the Walker breakdown are mandatory to draw conclusions about the 
phenomenological Gilbert damping parameter tuning the magnetisation dynamics. 

PACS numbers: 33.55.Fi, 75.70.-i, 42.65.-k, 75.30.Gw 
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I. INTRODUCTION 

Analysing the behaviour of a magnetic domain wall 
under the action of an external magnetic field is not a 
new concept. For example, magnetic bubble memories, in 
which information was stored in small magnetised areas 
delimited by domain walls and written by displacement 
of the walls under magnetic fields, were designed almost 
30 years ago Recently, an alternative approach has 
been proposed for reversing the magnetisation of a mag- 
netic cell, namely current-induced magnetisation rever- 
sal 0, In this frame, the use of current-induced mag- 
netic domain wall propagation, as suggested by the pio- 
neering work of Berger [3,, Jj] is nowadays raising a grow- 
ing interest. Controlling domain wall motion induced by 
a magnetic field and/or by a current flowing directly into 
the magnetic element, is thus a main issue in magnetism 
as well as a technological challenge. Indeed, in the at- 
tempt of designing nanowires enabling fast domain wall 
motion, disperse results have beenjjublished so far in re- 
gards of the field-induced motion 
for the current-induced one [l^, [Dj, lU, I13|, 11 
Moreover, the tuning of wall damping and motion re- 
mains an open problem of magnetisation dynamics. 

This article intends to give a global vision of the de- 
pendence of domain wall velocity upon magnetic field and 
current density strengths in nanowires. We analyse the 
field and current induced propagation, within the frame- 
work of previous analytical studies for the viscous regime 
in thick infinite samples containing a single 180° Bloch 
domain wall driven by a magnetic field Pjllall^- Indeed, 
when the effect of the applied magnetic field is sufficient 
to overcome pinning forces, the wall velocity v depends 
linearly on the magnetic field H, corresponding to viscous 
motion. Two linear velocity regimes, separated by a com- 
plex transient regime, have been predicted to exist, below 




as well as 

il[il[ii 



and well above a threshold field called the Walker field 
Hw [3, delineating the stability of the wall structure. 
In each linear regime, a relation between the wall mobil- 
ity m = dv/dH, the intrinsic Gilbert damping parameter 
a [20| and th e dy namic domain wall width has been es- 
tablished [H, Ef . There is in the literature a significant 
dispersion in the interpretation of velocity measurements 
in nanowires [l,0j[l,li3| i particularly in relation to damp- 
ing parameters extracted from the latter. This spread, 
partly addressed in [8] , shows how crucial the knowledge 
of the Walker field is. We demonstrate that the Walker 
field, calculated for the case of infinite and thick samples 
containing Bloch walls, is reduced in the case of confined 
systems. We find also that the same geometry factor af- 
fects the Walker-like critical current density which limits 
the stability regime of the domain wall under current in 
nanowires. 

Up to now, the most widely used magnetic layers and 
devices rely on domain wall motion in systems with an in- 
plane in^netic anisotropy d, i, 0, i, H 0, EH, El, El, E3, 
El, El, E3 ■ However media with out-of-plane anisotropy 
are promising candidates for improving the storage den- 
sity and for designing magnetic logic devices f2ll| . In this 
work, perpendicularly magnetised systems will be exten- 
sively investigated; however the presence of an in-plane 
anisotropy will not affect the geometrical parameters rel- 
evant for the wall motion outside of the pinning limited 
regime. 



II. CURRENT AND FIELD INDUCED 
PROPAGATION OF A BLOCH WALL 

We focus on the behaviour of the local magnetisation 
M inside a domain wall, oriented as defined in Fig. ([T^) . 
The length of the out-of-plane magnetised track sketched 
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in Fig. is supposed to be infinite, and oriented along 
y. One first tries to evaluate the velocity of a wall prop- 
agating towards the positive y values under the action of 
a magnetic field H applied along the easy axis of mag- 
netisation z and/or a current density Je along y. We 
focus on structures with a 180° Bloch domain wall. 9{y) 
given in Eq. (fTa|) stands as a possible wall profile; the wall 
center is located at yo where = ^ and the characteris- 
tic domain wall width A is given in Eq. (jlbp . Although 
magnetostatic effects within domain walls are naturally 
taken into account in numerical work, they are rarely 
discussed in analytic treatments [l^]. In Eq. (jlbp and 
in the following, the demagnetising factors N^, Ny and 
Nz relate to the domain wall itself, Ms is the saturation 
magnetisation, A and Ki are the exchange and intrinsic 
anisotropy constants (see methods). 




b) Perpendicular anisotropy 
Z 



6{y) = 2 arctan 
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A 



K, + 2t:M^ {N^ cos2 (p + Ny sin^ (p - N^) 



(la) 
(lb) 



For (y5 = or TT, one has an "ideal" Bloch wall with 
A„=n = / A This expression involves 

the wall demagnetising factors along those directions in 
which the system has reduced dimensions: along the 
wire width and through the layer thickness (see meth- 
ods). 

Besides the shape and size of the sample, among the 
parameters affecting its mobility, the domain wall struc- 
ture is often neglected in analytic calculations. In partic- 
ular, for in-plane anisotropy, several wall types may exist 
depending on the film geometry: C-shaped or Neel like 
transverse walls in infinite films [1^ ; Neel-like transverse 
walls or vortex walls for a narrow wire ^23|. The system 
of coordinates sketched on Fig. ([TJ;) allows one to adapt, 
via a permutation of the coordinates, the calculations 
made here for Bloch walls in out-of-plane spin configura- 
tions to transverse walls in in-plane spin configurations. 
Vortex walls, frequently observed in experiments on large 
nanowires, are not considered in the present article, nor 
edge effects. 

Magnetisation dynamics are described by Eq. in 
which the precession of the magnetisation proposed by 
Landau and Lifshitz [2^, the dissipative effects intro- 
duced by Gilbert and the usually admitted adia- 
batic and non-adiabatic spin torque terms [1^, [13, [l^ 
reflecting the action of current, assuming the geometry 
of Fig. ((TJj) are merged together. Discussion on micro- 
scopic spin polarised electronic transport is beyond the 
scope of this paper, in which we just discuss the macro- 
scopic co nseq uences of the presence of the spin-torque 
terms [H,!!!!!!. 



c) Planar anisotropy 



FIG. 1: a): Polar d and azimuthal </p angles and the associ- 
ated spherical coordinate system (er,eg,e^) defining the ori- 
entation of the magnetisation M relative to a Cartesian axis 
system, b) Sketch of a 180° domain wall in a track of width w 
and thickness t showing a perpendicular magnetic anisotropy. 
c) Sketch of the system of coordinates allowing to adapt the 
calculations to a narrow track with a planar anisotropy. A 
denotes the wall width and H the magnetic field applied along 
z, defined as the easy magnetisation axis of the magnetic layer 
in both cases. The current density Je along y is considered as 
positive if the electrons flow towards the positive y values. 



dM 



dt 



= jHeff X AI- 



a dM dM 



f3u - dM 
(2) 

In this expression, 7 is the gyromagnetic ratio (here a 
positive quantity) and i?e// is the total effective field. P 
is the phenomenological non-adiabatic spin transfer pa- 
rameter 26, 27, 28]. u (Eq. [3]) has the dimension of a 
velocity and scales as the electrical current density Je'. 



2eM, 



(3) 



g is the Lande factor, /is is the Bohr magnetron, e is the 
electron charge and P is the polarisation factor of the 
current. 

We express the different torques acting on the mag- 
netisation of the wall region in the (eV,ee,e(p) spherical 
system of coordinates (see methods). From the general 
kinetic momentum theorem, it is possible to relate the 
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torque components Tg and and their associated pre- 
cession velocities 9 and ip: 



III. WALL VELOCITY AND STABILITY IN 
SYSTEMS OF REDUCED DIMENSIONS 



dt 

dip 



Ms 



(4a) 
(4b) 



The first torque to be considered (see methods) is that 
of the external field namely Vu ~ M x H (Eq. [5|). 
The second one Th^ = M x Hd (Eq. [5]) is that of the 
demagnetising field within the wall. The third torque 
r^r^ = M X Ha (Eq. [7]) is that of the equivalent damping 

field 



For a low applied field or current density, in ultra- 
thin perpendicular magnetised systems, pinning domi- 
nates and domain walls move following the creep law [2^ . 
In this non-linear regime, the motion is thermally acti- 
vated and proceeds by discrete jumps from one pinned 
configuration to the next. In the following, we assume 
that pinning effects are overcome and consider only vis- 
cous propagation. We first consider the steady domain 
wall motion, whose signature is a wall moving with a 
time independent azimuthal angle ip: ip = which im- 
plies that Tip (Eq. I10b|l must vanish. Writing = for 
6 — Tr/2 (at the wall center where Fi^ has an extremum) 
gives: 



H 






' Ms H sine 



(5) 



= 4ttM II (Ny - N^) sin sin v? cos | (6) 

sin 9 cos 9[Nz — Ny sin^ ip — cos^ ip] 




(7) 



In the same spirit, assuming that the profile at rest 9(y) 
given in Eq. [Ta|is conserved under field and/or current. 



which gives 



dy 



the adiabatic F^ and the non- 



adiabatic F^ spin transfer torques are expressed by Eq.[8] 
and El 



udMd9 M,usm9_ 



F« = 



7 d9 dy 



7 A 
dM 89 



(8) 



The total polar and azimuthal torques Tg and F;^ are 
given in Eq. llOal and Eq. llObI respectively. 



Tg = A'KMj{Ny - N^) sin 9 sin ip cos ip (10a) 



aMs . . ^ MsU sin0 
ipsm 6 H — 



sm2(p= ; ■ — - (11) 

only valid if \H + {(3 ~ a);^| < Hw (12) 

where we define the adapted Walker field by 

Hw = 2naMs\Ny - N^\ (13) 

Eq. [Insets a limit on the combined strength of the ap- 
plied field and current density for steady motion to occur. 



Hp 



i/3-a) 



■yA 



determines the Walker-like stability 

condition for the wall motion as a whole. This limit con- 
dition is strictly equivalent to the Walker breakdown con- 
dition, in the case where only an external magnetic field 
is applied. Above the Walker breakdown, the stability 
condition = is broken and the motion is no longer 
steady. There, the torque of the external field/current 
dominates those of the demagnetising (Eq. [S]) and damp- 
ing (Eq. [7|) fields. It is the ip torque component which 
fixes the steady regime limit for a wall but the under- 
lying quantity governing the wall's stability is the mag- 
netostatic field inside the wall, making its motion geom- 
etry dependant. As a direct illustration, in the steady 
regime (ip—O) and without current, the only non zero 
component of the damping torque is the ip component, 
which is proportional to the rate of change of 9 and de- 
pends only on the demagnetising field. As com par ed to 
the expression obtained by Schryer and Walker [l^, our 
expression of Hw (Eq. [T31) includes a geometry factor 
I Ny — Nx I . For a zero applied field {H = 0) one similarly 
gets a Walker breakdown current density Jw (Eq. [H)) . 
The same geometrical reduction factor enters the critical 
Walker-like current density defining the stability condi- 
tion of the wall. 



r<p = -MsH sin 9 



aMs 



MgU sin 9 



7 7 A 

AnMs sin 9 cos 9[N^ - Ny sin^ ip - cos^ ip] 



(10b) 



Jw 



AiraeM^ 7A 



\Ny - 



(14) 



The wall velocity is related to the way the spin aligns 
along the external field, i.e. to the rate of change of 9 and 
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therefore comes out from the resulting torque Tg given in 
Eq. llOal Making use of Eq. [H in the steady regime {(p = 
0), 9 is simply —7 x AnMsiNy — N^) sin 9 simp cos (p — 
ysme^ The velocity v of any spin embedded in a Bloch 
wah is related to the characteristic width by: v — — ^9. 

sin & 

Choosing again 9 — tt/2 and using Eq. [Til we obtain: 



Walker breakdown 



7A 

'^steady — 



H+iP-a) — 



7A 



H 



(15) 

This expression is formally equivalent to the one ob- 
tained in the Walker theory the magnetic field be- 
ing here shifted by a spin transfer equivalent field (3u/jA. 
This extra term is linked to the non-adiabatic spin trans- 
fer, which is the only spin transfer component efficient for 
steady motion. As shown for in-plane magnetised sam- 
ples, the adiabatic component distorts the wall structure 
but does not induce a stationary motion [H, [2I] . It is 
worth noting that the "/3 term" does not appear explicitly 
in the Tg torque component. It only appears implicitly 
via its action on the (p precession. 

li \H + {(3 - a)::^\ > 2TTaMs\Ny - N^], we assume 
that the precession of the magnetisation occurs at a con- 
stant angular speed so that the average value over time of 
sin2(p{yo, t) can be set to zero. Thus, averaging Eq. llOal 
and Eg. II Obi over a period of the precession of ip and mak- 
ing use of Eq. IH one gets the average velocity far above 
the Walker breakdown, simply given by Eq. 1161 



7A- 



H 



(3u 
7A 



(16) 



This regime is similar to the usual high field one de- 
scribed for a 180° Bloch wall [T, Tsl, fl3|, i.e. the av- 
erage velocity is linear in the field, following an initial 
drop in the mobility at the breakdown. The current ac- 
tion is thus i) equivalent to an additional magnetic field 
f3u/jA applied along the anisotropy direction of the sys- 
tem [2^, and ii) an extra velocity linked to the adi- 
abatic spin transfer, that is u/{l + a^), far above the 
Walker breakdown. Both linear velocity regimes are de- 
picted in Fig. [21 under sole action of a driving field (a) 
and under the combined action of a field and a current 
(b). 

In a one dimensional statement of wall motion, with- 
out current, a qualitative understanding of the damp- 
ing/demagnetising interrelation is easy to provide. When 
a drive field is applied along the anisotropy direction, 
M starts a ip precession movement, and tilts away with 
respect to its equilibrium orientation at rest, according 
to the external field torque given in Eq. [5] The out of 
the xz wall plane component of the magnetisation cre- 
ates magnetic charges. The 9 component of the torque 
resulting from the induced magnetostatic field (Eq. [S]) 
describes the resulting additional precession of the spin 
around Hd- Damping and demagnetising torques then 
rapidly "freeze" the precession by cancelling the Zeeman 




a) 



precessional 
propagation 



Under field 



27iaMs|Ny-Nx| 



H 




b) 



/(l+a) 



Under field and current 



-au/yA iTiaMs |Ny - Nx| H +((3-a) u/yA 

FIG. 2: Sketch of a 180° domain wall's velocity as a func- 
tion of (a) an external magnetic field H and (b) a driving 
term combining a field and a current {u relates to the current 
density) in a nanowire. This cartoon indicates the two linear 
regimes of velocity, below and far above the Walker break- 
down. The dotted line in the transient non-linear regime is a 
guide for the eyes. 



torque. At the Walker field, the equilibrium that ex- 
ists between the damping/demagnetising torque and the 
drive field torque keeping ip constant is broken and the 
wall exhibits its maximum velocity for = 7r/4 in the lin- 
ear region. Above the Walker breakdown, the azimuthal 
angle precesses around the external field. The motion is 
oscillatory, due to the periodic nature of the torque term 
with respect to p. In the case where there would be no 
damping, the wall would move back and forth about its 
initial position. In the presence of damping, the forward 
displacement is larger than the backward one. This re- 
sults in a positive average displacement during one period 
and hence in a non zero average velocity. 



IV. DISCUSSION 

In contrast to the case of a thick infinite film, one 
must consider all components of the wall's internal de- 
magnetising field to describe the Walker breakdown of 
the wall under the action of either a magnetic field or 
a spin polarised electrical current when the system's di- 
mensions are comparable to, or less than, the domain wall 
width. For a thick and infinite sample, the only non-zero 
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wall demagnetising coefficient is the one along y since 
the magnetic charges along z (thickness) and x (width) 
are repelled towards infinity (Fig. [3]) . Said another way, 
in that geometry, the demagnetising field has only one 
significant component which is perpendicular to the wall 
plane and the additional torque related to Hd^ (through 
the wall width) is the only one that is efficient in align- 
ing the spin and pushing the wall forward. For a track 
of reduced dimensions and/or an ultrathin film (t -C A, 
ui <C A), the coefficients N.j. and are not zero any- 
more and reflect the aspect ratio w/A and t/A (Fig. [3]). 
Since + Ny + Nz — 1, the demagnetising field Hdy 
perpendicular to the wall and its associated torque will 
be reduced with respect to the infinite thick film case. 
Time resolved Kerr microscopy experiments performed 
on continuous Pt (4.5 nm) / Co (0.4 nm < i < 0.9 nm) / 
Pt (3.5 nm) layers grown by sputtering with a perpen- 
dicular to the plane magnetic anisotropy support the 
present conclusions [3Q] . In the case of a nanowire of 
width u'=200 nm patterned from such a Pt(4.5nm)/Co 
(0.5 nm)/Pt(3.5nm) film, we model the ID wall cross- 
section as an ellipse which allows use of the standard 
expression Ny ^t/{t + A) and ^t/{t + w) [3l|. This 
leads to ~ 0.0025, Ny ~ 0.0775 and N, ~ 0.92. The 
Walker threshold field Hw of the nanowire is then about 
7.5% of that of the corresponding thick and infinite layer. 



a) Thick fllm 



A (« t) 




b) Thin fllm 



A(>0 



FIG. 3: Sketch of tlie yz plane of a 180° domain wall with 
magnetic charges and the associated y component of the de- 
magnetising field (Hdy) through the wall plane for different 
thicknesses. Hd is much smaller in case b) than in case a). 



plane magnetisation) of thickness ranging between 310 
and 3000 A [51]. The authors verified the Bloch nature 
of the domain wall (above 1000 A) and then observed an 
increase of the mobility with the thickness accounted by 
the domain wall width evolution; for smaller thicknesses 
and/or samples containing Neel walls, either it was not 
possible to observe the viscous regime or the thickness 
dependance was inverse. In the case of a few hundreds 
A thick NiFe nanowires of sub-micrometric width, v{H) 
experiments have been reported d, 0, H, However, 
the two successive regimes and the Walker threshold have 
been evidenced only recently in in-plane systems [1, [lo| . 
For current induced motion, the Walker breakdown was 
also evidenced recently Satisfactorily, in [ll[, 

the damping constant used to account for the experimen- 
tal value of the Walker field is compatible with our state- 
ment about Hw. 



V. CONCLUSION 

In summary, in the frame of designing devices based 
on domain wall motion, understanding of the Walker 
breakdown is critical. Above this breakdown, the wall 
structure is unstable causing the instantaneous wall 
velocity to oscillate and the average velocity to drop 
abruptly. We have demonstrated that wall confinement 
effects must be taken into account. In nanostructures, 
magnetostatic effects within the wall lead to a restriction 
of the stability region relative to thick and infinite films 
(Eq. [T5)) ; the domain wall width parameter A (Eq. Ilbp 
is also modified. The domain wall structure tuning 
is also a tricky problem and no consensus has been 
achieved yet concerning the evaluation of a. For current 
induced domain wall motion, dispersed results have 
been obtained so far, partly because of the distortion of 
the wall, especially in in-plane systems. Therefore, the 
choice of perpendicularly magnetised ultrathin layers, 
whose wall structure is known and in which the non 
adiabatic contribution should be large [l^ , are probably 
good candidates for a detailed understanding of the 
underlying mechanisms of wall motion and damping as 
well as for reliability of devices. 



It is in principle possible to evaluate the damping con- 
stant a from domain wall velocity experiments. Under 
field, this has been tentatively done using in-plane mag- 



netised samples 
magnetised systems [32 




but rarely for out-of-plane 
33|, in which the creep regime 
often extends over a large field range [1^ and prevents 
the observation of the linear regimes. The velocity de- 
pendence over an extended field/current range and the 
knowledge of the Walker breakdown are mandatory to 
draw conclusions about the phenomenological Gilbert 
damping parameter a. Experimentally, the first unam- 
biguous evidence of the initial linear regime has been 
obtained in the case of large area NiFe films (with in- 



VI. METHODS 
A. Wall's parameters 

In the case of a uniformly magnetised ellipsoid, the in- 
ternal magnetostatic field Hd is proportional to the mag- 

/ \ 
netisation according to Hd = — 47r Ny M. 

V N, J 
The demagnetising factors N^, Ny and Nz along the x, 
y and z directions respectively verify N^ + Ny + N^ = 1. 
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Magnetostatic demagnetising factors are complex to cal- 
culate in general l34l| but they can be evaluated for 
simple shapes [3ll.l34{. 

In the case where an out-of-plane uniform magnetic 
layer (and not a domain wall) is considered, an effective 
anisotropy constant is measured: K — Ki — 27rAf^, that 
takes into account the intrinsic anisotropy constant Ki 
and the demagnetising field of the out-of-plane saturated 
layer [Nl^y'''' = N^^v"'' = and Nif'^''^ 1). 
In the case of a pure Bloch wall {(f = 0), one can com- 
pare A^^o with the domain wall width parameter of a 
Bloch wall at rest in an infinite and thick layer p^ . 
that is ^*'Bioch = A/Ki. For an infinite but thin 
layer, the Bloch domain wall width becomes I^biocK — 
^A/{Ki- 27rM^N,) = ^A/iK + 2^M]Ny). 



B. Torque calculations 

The torques of the exchange and anisotropy fields com- 
pensate each other and therefore are not included in the 



torque calculations. They define the wall profile at rest, 
that is considered as conserved when investigating the 
dynamics. 

The active torques given in Eq. [5l [5] and [7] 

are obtained using the transfer matrix Yr,e,ip = 
sin 6 cos if sin 6 sin if cos 9 
cos 6 cos if cos 9 sin ip — sin ( 
— sin if cos if 
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